Creep failures in heterogeneous materials 
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We present creep experiments on fiber composite materials with controlled heterogeneity. Recorded strain 
rates and acoustic emission rates exhibit a power law relaxation in the primary creep regime (Andrade law) 
followed by a power law acceleration up to rupture over up to four decades in time. We discover that the failure 
time is proportional to the duration of the primary creep regime, showing the interplay between the two regimes 
and offering a method of rupture prediction. These experimental results are rationalized by a mean-field model 
of representative elements with nonlinear visco-elastic rheology and with a large heterogeneity of strengths. 
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The damage and fracture of materials are technolog- 
ically of enormous interest due to their economic and 
human cost. Failure of composite systems is particu- 
larly important in naval, aeronautics and space indus- 
try. Despite considerable experimental {]]> an d 
theoretical work |4, Sla] on fracture, many questions 
have not been answered yet. Recently, statistical physi- 
cists have shown the existence of a power law acceler- 
ation of acoustic emissions announcing the global fail- 
ure of heterogeneous materials 0, Hf, similar to the 
critical behavior of out-of-equilibrium phase transition 
|6], offering a way to predict material failure |2]. 

This letter presents creep experiments on compos- 
ite materials, which we explain using a simple model 
of representative elements, in the framework of fiber 
bundles models. Creep is the progressive deformation 
of a material under constant load at a given tempera- 
ture. Three creep regimes are usually observed. During 
the primary regime, the strain rate decays as a power 
law with time following the application of the stress 
(Andrade law) |7]. The secondary regime describes 
a quasi-constant deformation rate, which evolves to- 
wards the tertiary creep regime, if the stress and the 
temperature are high enough, during which the strain 
rate accelerates up to rupture. 

The experiments are carried out on cross ply 
glass/polyester composite materials and on Sheet 
Molding Compound (SMC) composites. Two type 
of cross angle ply laminates are fabricated, denoted 
[±62°] and [90°/35°], where the angles measure the 
directions of the glass fibers with respect to the load- 
ing direction, with a fibre volume fraction of 75%. The 
SMC composites consist of a combination of polyester 
resin, calcium carbonate filler, thermoplastic additive 
and random oriented short glass fibers, in the form of a 



sheet. The relatively low fibre volume fraction, about 
30%, and the uncontrolled filler and reinforcement dis- 
tribution during processing lead to a more heteroge- 
neous structure for the SMC compared to the cross ply 
composites. The [±62°] and [90°/35°] specimens have 
dimensions 14 x 100 x 2 mm 3 . The SMC samples are 
in the form of 120 mm barbell with 3 mm tickness. All 
specimens are subjected to a constant stress s and tem- 
perature T (below the glass transition of the matrix), 
which were fixed to s = 15 MPa and T = 60°C for 
the [±62°] specimens, s = 22 MPa and T = 60°C 
for the [90°/35°] specimens, and s — 48 MPa and 
T = 100°C for the SMC. The creep tensile tests were 
performed using a servo-hydraulic mechanical testing 
system. Constant tensile load was applied and the re- 
sulting strain and acoustic emissions were recorded un- 
til final rupture. Acoustic emissions (AE) is a standard 
technique to monitor the evolution of damage in com- 
posites, due to matrix cracks, fiber matrix debounding, 
fiber breaks, and delaminations |9]. We used a Mistras 
data acquisition system by Physical Acoustics Corpo- 
ration with 2 resonant sensors (200 kHz - 1 MHz). 

Normal primary creep transients followed by sec- 
ondary and tertiary creep were observed for almost all 
samples, both for the strain rate (Fig. ^ an d f° r the 
AE rate (Fig. |2j. The decrease of the strain rate and 
EA rate in the primary creep regime can be described 
by Andrade's law |7] de/dt ~ t~ p , with an exponent p 
in the range 0.2-1.4 for the 15 samples tested |8]. The 
crossover for small times is probably due to the fact that 
the stress progressively increases up to about 10 sec af- 
ter the start of the experiment. A quasi-constant strain 
rate (steady-state or secondary creep) is observed over 
an important part of the total creep time, followed by 
an increase of the creep rate up to failure in the tertiary 
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creep regime. Creep strains at fracture are around 40% 
for angle cross ply composites and « 4% for the SMC. 
The acceleration of the strain rate before failure is well 
fitted by a power-law singularity de/dt ~ (t c — t)~ p 
with p' in the range 0.3-1.1 depending on the sample 
|H • The critical time t c determined from the fit of the 
data with a power-law is close to the observed failure 
time. Our experiments confirm over large time scales 
covering up to four orders of magnitude in time previ- 
ous announcement of power laws in the tertiary creep 
regime, which were established over more limited time 
scales |3]. We also obtain the same temporal evolu- 
tion for the AE energy rate, with larger fluctuations for 
the energy rate than for the event rate due to the exis- 
tence of a power-law distribution of AE energies. We 
found no qualitative differences between the behavior 
of the [±62°] and [90°/35°] cross ply composites. The 
values of p and p 1 are on average a little larger for the 
SMC than for the cross ply composites, possibly due to 
the larger heterogeneity of the SMC or to the different 
values of the applied stress. 

There is a huge variability of the failure time from 
one sample to another one, for the same applied stress, 
as shown in Fig. [3] This figure shows that the transi- 
tion time t m between the primary creep regime and the 
tertiary regime, measured by the minimum of the strain 
rate, is proportional to the rupture time t m as 2/3 t r . 
We also found a negative correlation between the An- 
drade exponent p and the rupture time t r 0]. These 
observations shows that damage in the primary regime 
impacts on its subsequence evolution in the secondary 
and tertiary regimes. This suggests a way to predict 
the failure time from the observation of the strain rate 
or AE rate during the primary and secondary creep 
regimes, before the acceleration of the damage leading 
to rupture. 

Creep observations have been modelled in terms of 
visco-elastic fibers, with deterministic dynamics and 
quenched disorder 1 1 ] . This model reproduces a 
power law singularity of the strain rate before fail- 
ure with p' = 1/2 in the case of a uniform distribu- 
tion of strengths 1 1 0] but does not explain Andrade's 
law for the primary creep. Here, we start from the 
model of [10] and enrich it with a more realistic rhe- 
ology and heterogeneity, in order to account simulta- 
neously for Andrade's law in the primary creep and 
for the power-law singularity of the strain rate before 
failure. We view a composite system as made of a 
large set of representative elements (RE), each element 
comprising many fibers with their interstitial matrix. 
The applied load is shared democratically between all 
RE. This assumption has been shown to be a good ap- 
proximation of the elastic load sharing for sufficiently 
heterogeneous materials [13]. Each RE is modelled 



as a non-linear Eyring dashpot lHUl in parallel with 
a linear spring of stiffness E. The Eyring rheology, 
which is standard for fiber composites, consists at the 
microscopic level in adapting to the matrix rheology 
the theory of reaction rates describing processes acti- 
vated by crossing potential barriers. A given RE fails 
when its elongation/deformation e reaches a thresh- 
old. The rupture thresholds are distributed according 
to the cumulative distribution P(e) given by P{e) — 
1 — ( e oi/( e + eo2))'\ where eoi and eo2 are two con- 
stants with eoi < eo2- The fraction 1 — (601/602)^ 
breaks as soon as the stress is applied. The power-law 
distribution P(e) for large e is motivated by the large 
distribution of failure times for the same applied stress 
(Fig. |3- The exponent /1 > 1 controls the amplitude 
of the frozen heterogeneity of the RE strengths. 

The equation controlling the deformation e(t) of 
each surviving RE is 

with the initial condition e(t = 0) = 0. The fraction of 
unbroken RE is 1 — P(e) and s/(l — P(e)) is the stress 
applied on each unbroken RE. 

The system defined by Q is stable (no global rup- 
ture) if the differential equation Q has a stationary 
solution de/dt = with e > 0, i.e., if the equation 
((e + eo2)/eoi) A ' = (E/s)e has a non-trivial solution. 
This defines a threshold s* below which the strain con- 
verges asymptotically to a constant and above which 
de/dt grows without up to rupture. 

In the primary regime e <C eo2 thus (e + eo2) M ~ 
eo2 M (l + A*e/eo2). If the stress on the dashpot is small, 
we can replace sinh by exp /2. With these approxima- 
tions, the differential equation Q has the solution 



* 2e-*(3Sr + tK(i(E-^(^Y) ' 

Expression (0 predicts that, if the stress is not too 
large, de/dt is of the Andrade form ~ t~' p , with an ex- 
ponent p — 1 at early times. For larger s, the strain rate 
starts to accelerate as soon as the load is applied. Note 
that the observation of Andrade's power-law creep in 
this model does not involve any failure of RE and is 
thus independent of the choice of the distribution of 
rupture thresholds P(e). 

In the tertiary creep regime, we can neglect eo2 com- 
pared with e. Close to failure, for large e, the lin- 
ear term Ee is small compared with s/(l — P(e)) — 
(e + eo2) M if H > 1. This leads to the equation 

de K />e"\ 
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Its solution is, to leading logarithmic order, 

where A = eoi (/3s) _1 . Contrary to the primary 
regime, the heterogeneity of the rupture threshold is an 
essential ingredient for the power-law singularity be- 
fore failure, but the leading power-law term with p' = 1 
in @ does not depend on the exponent fi characteriz- 
ing heterogeneity. The acceleration toward failure is 
due to the positive feedback effect of broken RE, which 
increases the stress and strain on the unbroken RE lead- 
ing to the global failure of the system. 

Figure |4] shows the numerical solution of equation 
Q together with the approximate analytical solutions 
in the primary creep and ® close to failure, for 
different values of the applied stress s. In the primary 
creep regime close to the rupture threshold s « s*, 
we observe numerically an apparent exponent p < 1, 
smaller than predicted by (0, which can explain the 
values of p found experimentally |8]. For a stress 
s ^ s*, the strain rate accelerates immediately when 
the load is applied. For s < s*, the p- value decreases 
in the model between 1 and as the applied stress in- 
creases. The duration of the primary creep also de- 
creases with s. The model thus explains the correla- 
tion found experimentally between the p-value and the 
failure time 1 8 ] . 

In the tertiary regime, for s 3> s*, we find numeri- 
cally that expression Q is a good approximation very 
close to failure t « t c . But for s sa s*, there is a 
crossover further from failure with an apparent expo- 
nent p' = 0.9. This simple model thus reproduces both 
power-laws in the primary and tertiary creep regimes, 
with an apparent exponent p < 1 for the primary creep, 
and with p' = 1 for the tertiary regime, except for a 
crossover with an apparent exponent p' a little smaller 
than 1 . This crossover with p' < 1 is however not suf- 
ficient to explain the observations of p' = 0.93 over 4 
orders of magnitude in time t c — t (Fig. Q. 

The failure time has a power-law singularity ~ (s — 
s *-j-i/2 f or s ^ s *, as found previously 11011 for the 
model with a linear dashpot, and decays exponentially 
for s ^> s* 1 8]. The transition time t m (minima of the 
strain rate) is equal to t c /2. This result recovers the 
proportionality of t m and t c found experimentally, but 
predicts a duration for the primary creep shorter than 
the observations t m « 2i c /3 (Fig.|5|l. 

In conclusion, we have shown that the interactions 
between the RE elements together with a large hetero- 
geneity and a simple nonlinear rheology is sufficient to 
explain qualitatively and quasi-quantitatively our ex- 
periments. This model replaces the need for com- 
plex memory effects (such as the integro-differential 



Schapery long-memory formalism 1 12]) often invoked 
in the composite literature. A natural improvement of 
the model would be to relax the democratic load shar- 
ing rule as in 11011 in order to introduce realistic elastic 
interactions. This improvement may provide a more re- 
alistic value of p' and of the constant of proportionality 
between t m and t c (Fig. [3}. 
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FIG. 1: Creep strain rate for a [±62°] specimen, (a) linear time scale, (b) logarithmic time scale to test for the existence of 
Andrade's law in the primary creep, (c) logarithmic time scale in t c — t to test the time-to-failure power law in the tertiary creep. 
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FIG. 3: Relation between the time t m of the minima of the strain rate and the rupture time t r , for all samples. 
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FIG. 4: Strain rate de/dt given by Q for different values of the stress s, and with parameters E = 20 GPa, fi = 1.2, 
ecu = 0.003, eca = 0.015, /3 = 50 GPa -1 and K = 10~ J sec - . Panel (b) illustrates Andrade's law in the primary regime, 
with exponent p ~ 1 for s = 22 MPa and p ~ 0.8 for s = 25 MPa. The dashed line is the approximate solution J2j of Q 
with s = 22 MPa. Panel (c) shows the power law acceleration of de/dt before failure for s — 25 MPa and s — 33 MPa, with 
p' ~ 1 asymptotically. In (a) and (c) the time is normalized by the rupture time for s = 25 MPa and for s = 33 MPa, and by 
the time when de/dt decreases below 10 -14 for s — 22 MPa. 



